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Abstract 

Interference alignment (IA), given uncorrelated channel components and perfect channel state 
information, obtains the maximum degrees of freedom in an interference channel. Little is known, 
however, about how the sum rate of IA behaves at finite transmit power, with imperfect channel 
state information, or antenna correlation. This paper provides an approximate closed-form signal- 
to-interference-plus-noise-ratio (SINR) expression for IA over multiple-input-multiple-output (MIMO) 
channels with imperfect channel state information and transmit antenna correlation. Assuming linear 
processing at the transmitters and zero-forcing receivers, random matrix theory tools are utilized to 
derive an approximation for the post-processing SINR distribution of each stream for each user. Perfect 
channel knowledge and i.i.d. channel coefficients constitute special cases. This SINR distribution not 
only allows easy calculation of useful performance metrics like sum rate and symbol error rate, but 
also permits a realistic comparison of IA with other transmission techniques. More specifically, IA is 
compared with spatial multiplexing and beamforming and it is shown that IA may not be optimal for 
some performance criteria. 

I. Introduction 

Many important wireless communication scenarios can be modeled using an interference 
channel. Examples include interfering base stations in a cellular network, wireless local area 
networks, and simultaneous transmission in mobile ad-hoc networks. The general capacity of 
the interference channel and the design of practical schemes approaching the known upper 
bounds on sum rates have been of great interest over the last 30 years. The earliest attempts to 
characterize the capacity region of the interference channel, inspired by the framework established 
by Shannon in [0, were focused on two-user interference channels. Although the special cases 
of strong and very strong interference have been solved (3), flU, the general capacity of the 
interference channel is still an open problem. Recently, a series of attempts have been made 
to describe an approximation of the asymptotic sum capacity behavior known as the maximum 
achievable multiplexing gain or degrees of freedom (DoF) [5] where focus is on the high SNR 
regime and the interference/broadcast characteristics of the wireless network. The DoF studies 
paved the way for a novel method of dealing with interference, known as interference alignment 

(ia) a m. 



A. Recent Work and Motivation 

IA uses beamforming matrices at the transmitters to align the interference to a received signal 
subspace such that an interference free subspace becomes available for direct signal transmission. 
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Thus in IA, the primary goal of the transmitting nodes is to reduce their decremental effect 
on the unintended receivers. In contrast to other techniques of interference management such 
as orthogonal access, decoding the interference or treating the interference as noise 10, 
IA achieves the maximum DoF in a K-user interference channel [5|. IA has also been used 
successfully to characterize the maximum DoF of other network scenarios including the MIMO 
interference channel with time varying/frequency selective channels flU or constant channel 
coefficients ll9l iflOl the X channel Qj], the MIMO X channel flU, X Q2, and MIMO X 
networks [fT3l . 

Given the potential of IA, recent work has further explored its applications and limitations. 
Iterative and distributed algorithms for IA over constant channel coefficients are presented in 
||9l , lfl4l . The feasibility of IA over spatial dimensions is studied in ifTOl . [fT5l . Reducing the 
overhead associated with IA is considered in ifToll . Adaptation of IA to multi-cell networks is 
considered in [fTTl , |fT8l . A method of opportunistic access in cognitive radios inspired from IA is 
proposed in |fT9ll . Using IA for secure communications is studied in GUI . And finally, extending 
IA to relay-aided networks is considered in [|2Ti - ll23ll . 

On the one hand, despite the large number of practical wireless networks for which IA is 
being considered, IA in its original form Q is only shown to be optimum for asymptotically 
high SNR given perfect channel state information (CSI) and i.i.d. channel coefficient values 
among the users and between the antennas. On the other hand, in practical systems E4l . the 
communicating nodes only have access to an imperfect estimate of the channel coefficients, 
working at intermediate SNR values is inevitable and correlation either between the nodes or 
especially between the antennas exists. Previous work to quantify the effect of imperfect CSI and 
intermediate SNR values on IA are few, and either confined to calculating bounds on achievable 
sum rates [fTTl or rely on experimental results ||23l without providing accurate quantification of 
key performance measures such as symbol error rate or achievable sum rate. The reason is the 
complicated expressions of the linear beamforming/combining filters. 

Our motivation for considering ZF receivers is two fold. First, although zero-forcing (ZF) 
receivers are asymptotically sufficient for achieving the DoF promised by IA [5 ] and they provide 
a simple and effective method for multiple stream detection, in the context of IA, little is known 
about the performance of such receivers in intermediate SNR regimes with imperfect CSI and 
correlated antennas/channels. Second, some other classes of receiver filters, such as minimum 
mean square error or regularized ZF filters, converge to the ZF receiver at high SNR which 
raises further question about ZF receivers at lower SNR. 

Note that interference can be aligned on both the signal level space ll25l (multi-level and 
lattice coding) and in the signal vector space 0, [fT3l . Il26ll (designing the beamforming and 
combining filters using the time/frequency/spatial dimensions). Furthermore, although next gen- 
eration MIMO based networks [24| present the opportunity to code over frequency, time and 
space dimensions, lack of sensitive frequency synchronization between the nodes, potential high 
correlation between the channel coefficient values (IA highly relies on independence between 
the used channel coefficients to achieve the full multiplexing gain [5]) and delay- sensitive 
communications justifies using IA only over the spatial dimensions, i.e. IA over constant MIMO 
channels ||9ll, IfTOll. 

B. Contributions 

In this paper, we quantify the IA performance under imperfect CSI for constant MIMO 
channels with transmit correlation with ZF filters. First, using random matrix theory tools, we 
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show that given perfect CSI at all the nodes and uncorrelated i.i.d. channel coefficients, the 
received SNR per stream for each user after a ZF equalizer is exponentially distributed. In 
other words, using a ZF receiver results in parallel single-input-single-output Rayleigh channels. 
Next, we show that if there exists an arbitrary Kronecker-modeled transmit correlation, using 
the asymptotic behavior of the eigenvectors of Wishart matrices, we can quantify its effect 
approximately on the received SNR distribution. Our analysis shows that the accuracy of this 
approximation depends on the number of antennas at each node, the transmit antenna correlation, 
and the transmit power. 

After quantifying the effect of transmit correlation on the SNR distribution, we analyze the 
effect of imperfect CSI (Gauss-Markov model) on the received SNR. We show that the imperfect 
channel knowledge reduces the mean received SNR in proportion to the number of streams. In 
other words, increasing the multiplexing gain, either by using more antennas at each node or by 
increasing the number of transmit-receive pairs, increases the detrimental impact of imperfect 
channel knowledge. Moreover, we show that if the imperfection in CSI does not vanish at 
asymptotically high transmit powers, one can not achieve the full multiplexing gain promised 
by IA. Finally we show that using the derived per-stream SNR distributions, it is possible to 
accurately compare an IA network with other interference management methodologies such as 
an orthogonal access network which is utilizing spatial multiplexing (SM) or beamforming. 
Our recent work on this subject HI did not include channel correlation and was restricted to 
quantifying the impact of imperfect CSI on the distribution of the post-processing SINR for IA 
over constant MIMO channels. 

C. Organization and Notation 

Organization: In Section [II] we present the system model and analyze the performance of IA 
given perfect channel knowledge and i.i.d. channel coefficient values. In Sections Hill and HVl we 
progressively introduce transmit correlation and imperfect CSI into the system model and then 
quantify their effects on the distribution of the post-processing SINR of each received stream. 
Section |V] analyzes a point-to-point MIMO system relative to our IA configuration. Numerical 
experiments are presented in Section [VI] and concluding remarks are given in Section IVII1 

Notation: Capital and small bold letters stand for matrices and vectors. A*, A T and A -1 are 
conjugate transpose, transpose and inverse of A respectively. A(n, m) is the element on the n th 
row and m th column of A. A(:,m) is the vector of m th column of A. tr(A) and rank(A) are 
trace and rank of A, respectively. [A, B] is the matrix constructed by horizontal concatenation 
of matrices A and B. vec (A) operator stacks the columns of A. \a] is the smallest integer 
larger than or equal to a, I N is the N x N identity matrix and axb is an a x b matrix of all 
zeros. 



Consider the f^-user N r x N t MIMO interference channel shown in Fig. \T\ Transmitter i 
encodes d{ streams using a precoding matrix F, which is then decoded by the i th receiver after 
processing the received signal with a combining matrix W». Assuming perfect synchronization 
between the nodes, the received signal at receiver i can be written as 



II. IA Over a Constant MIMO Channel 



A 




(1) 



k=l 
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where z, is AWGN with elements distributed as £/V(0, a 2 ), transmitted streams Xj obey the total 
power constraint E {x*Xj} = P, Hj^ is the matrix of the channel coefficients between receiver 
i and transmitter k. Note that in constant MIMO channels, coding over frequency/time does not 
change the achievable DoF [5] and without loss of generality we can assume the channels are 
narrowband. 

The goal of IA for the MIMO constant channel is to design the precoding and combining 
matrices such that the following conditions are satisfied 



rank(W i H ii F i ) = d; 
W,H ifc F fe = VA; ^ % 



\/ie{l,...,K}. 



(2) 



Except the 3-user [[51 and the (N + l)-user N x N [TOl constant MIMO channels, analytical 
expressions directly solving (O for the unknowns are still under investigation. There exist, 
however, numerical methods for both designing the precoding/combining matrices and calculating 
the achievable DoF [0, IfPfll . We focus in this paper on the alternating minimization method in 
||9l ; we leave extending the developed theory and the results to other IA precoding/combining 
matrix design techniques for future work. Note that a system of IA can have multiple solutions 
|fT5l and it is possible to select a solution with specific properties for each channel instance 
||27l Section VI]. Our analysis, however, is based on selecting a random IA solution for a given 
channel instance and studying the behavior of solutions with specific properties (through extreme 
statistics) is out of the scope of this work. 

For the i th transmitter/receiver pair, the alternating minimization method results in a unitary 
transmitter beamforming matrix, F i5 and a set of non-unique orthonormal basis for the interfer- 
ence subspace, i.e. columns of Cj, such that 



F*F, ; 



C*Q = I Vie {!,...,#}. 



(3) 



Note that in the alternating minimization method, the direct channel links Hjj do not appear 
in the computation of precoding matrices (and interference subspaces) flU and therefore the 
elements of Hjj are independent of Fj and Cj for % e {1, . . . , K}. 

The next step is to design the receiver equalizers. Assuming the alternating minimization 
method has converged to an IA solution (see [|9]| and [|28l Section 6] for discussions on its 
convergence), by knowing a basis for the interference subspace at each receiver i, i.e. columns 
of Q, ® can be written as 



K 



HjjFjX,; + Cj AjfcXfc + Zj = Hj 

k^i 



+ Zj 



(4) 



where A ik determines the interference from transmitter k at receiver i and is given by CjA ifc = 
H^Ffc and Hj = [HjjFj, Cj] is the effective channel at receiver i. Note that any other or- 
thonormal basis for the interference subspace is related to Cj through a unitary mapping and, 
by appropriately transforming A ik , we can restrict our attention to C, without affecting the 
forthcoming discussions. Moreover, as Jin is independent of Cj and Fj (through construction), 
the columns of the N r x di matrix of HjjFj do not completely lie in the subspace spanned by 
columns of the N r x (N r — di) matrix of Cj, i.e. the N r x N r matrix of [HjjFj, Cj] is full rank. 
In addition, given that IA is feasible and that transmitters are communicating the maximum 
number of allowed streams, Fj and Cj are always of dimension N t x di and N r x (iV r — di) and 
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Hj is always square. Therefore, the ZF equalizer at receiver i will be 



W? F = [I^O^-djHr 1 . (5) 

Define Bj = [i^, O^x^-di)] • Applying the ZF receiver given by © to the received signal 
given in © leads to an expression for the SNR of the n th stream at receiver i 



Hi,'. 



lo/di 



B,;fH?H,;)- 1 B? 



(6) 



n,n 



where 7 Q = Using Lemma [T} the denominator of © can be simplified into an expression 
better suited for statistical analysis. 

Lemma 1. Assuming IA is feasible, the denominator of© simplifies to [F*H* 4 (\n t — CjC*) HjjFj 
Proof: Given that Hj = [H^Fi, Q] 

[HjjFj, Cj] [HjjFj, Cj] = 



H*Hj 



r 3 



r 2 
r 4 



(7) 



where V, = F*H*jHjjFj, T 2 = F*IL* Q, F 3 = C*HjjFj and T 4 = C*Q = Noting that 

left and right multiplication of any matrix by Bj and B* keeps the first di rows and the first 
di columns of that matrix, if (T 1 — r 2 (r4) _1 T 3 ) is not singular, the inverse of © simplifies to 



(ri - r 2 (r 4 )- 1 r 3 ) 1 m Section 3.5.3]. As (I N , 
(Hermitian) idempotent matrix 

r p /T-i \ — 1-p I' ll II xp 1?*XJ* /""< f~ i *~f-J TT 

F*T_T* /'T f 1 I I 1/ 1?*XJ* p TJ 1? 



CjC*) is a projection matrix and hence an 



F - TT * p p TT TJ 



(8) 



where Cj = (l Nr — CjC*). Furthermore, Cj is a projection matrix into the di dimensional null- 
space of the interference subspace at the i th receiver constituting a combining matrix which 

di and the di x d$ matrix in © is nosingular. ■ 

lo/di 



satisfies ©■ Therefore, rank ^QHiiF; 
Using Lemma [H © simplifies to 

1i,n 



F*tj* r 1 1 1 1: 



(9) 



Define Aj = F^H^. Lemma [2] gives the distribution of AjA* which is then used in Lemma [3] 
to find the distribution of AjQA*. 

Lemma 2. If the channel coefficients are i.i.d. zero-mean unit-variance Gaussian random vari- 
ables, Hjj is independent o/Fj, (TJ]) holds, then the di x di random matrix AjA* is a complex 
Wishart matrix with N r degrees of freedom and covariance matrix 1^. 

Proof: In the alternating minimization algorithm, the columns of Fj are di eigenvectors (cor- 
responding to the di least significant eigenvalues) of the N t xN t matrix J^Lj ^ti (I ~~ CfcC*.) H ki . 

Let the eigenvalue decomposition of ^ Hj& (I — C*.C£) H*,j be UjSjU* where £j is a 
diagonal matrix holding the eigenvalues sorted in ascending order and Uj is a unitary matrix 
holding the corresponding eigenvectors. Then Fj = Uj [i^., O^x^t-*)] • Therefore Aj = 
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[1^, Odixim-di)] U?H£ and following the unitarily invariance property of H* i5 0^^^)] U*H*. 
has N r columns of size di with each element distributed as CjV(0, 1). It follows that the di x di 
matrix of AjA* has a central Wishart distribution with N r degrees of freedom and covariance 
matrix I d .. ■ 
Using Lemma [2] and noting that Q is a projection matrix into a sub-space of dimension di, 
Lemma [3] gives the distribution of AjQA*. 

Lemma 3. If the dxd matrix A A* has a central Wishart distribution with N degrees of freedom 
and covariance matrix 1^ and C is a projection matrix into a sub-space of dimension c < N, 
the dxd matrix of ACA* has a central Wishart distribution with c degrees of freedom and 
covariance matrix Id- 
Proof: As a projection matrix into a sub-space, C can be rewritten as USU*, where U 
is a unitary matrix and £ holds the c unity eigenvalues of C on its main diagonal and zeros 
elsewhere. Proof follows from unitarily invariance of columns of A and idempotent property of 
£ (S3 Theorem 3.4.4]. ■ 
From Lemma [3] it follows that the SNR of each stream given by © is exponentially distributed 
ED, i.e. 

/fa») = -exp(-^»Y (10) 



7o V 7o 

K 



By using a ZF receiver, the system effectively consists of \J d m parallel Rayleigh point-to-point 

m=l 

links. 



III. Transmit Antenna Correlation 



A more general channel model includes spatial correlation between the channel elements. In 
this section we suppose that there is transmit spatial correlation, thus 



H 



ik 



HW T> 



1/2 



\/i,ke{l,...,K}, 



(11) 



where LP ~ CjV(0, I) and R f is a constant Hermitian positive semidefinite (PSD) matrix 
11321 . Note that when the receive spatial correlation, R r , is not an identity matrix, neither the 
columns nor the rows of H& are independent and analyzing the distribution of H^CjHjj (which 
appears in the denominator of SNR expressions) is more challenging. Therefore, including receive 
correlation is a topic of future work. For the initial analysis, we assume that the spatial correlation 
is the same for all channels in the interference network; we relax this assumption at the end of 
the section. 

Assume IA is done over the instantaneous channels as before. Note that the feasibility of a 
system of IA is a function of number of equations and variables in © [fT5l . Therefore, as long as 
Rf is not rank deficient, transmit correlation will not change the achievability of IA. Replacing 
(flTT) in © gives 

lo/di jo/di 



F*(R, 1/2 )*(H-)*QH-RrF 



1/2, 



AjQA" 



(12) 
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~ 1/2 " ~ 

where Aj = F*(R/ )*(H^)*. A, has i.i.d. columns with covariance matrix Rj defined as 

R = covCAiC:, n), £,(:, n)) = E {f*(R 4 1/2 )* (H»(:, n))*H^(:, n)R] /2 F t ] = E {F*K t F t } . 

(13) 

Assuming Rj is known, similar to (PTOl) . the SINR of the n th stream at receiver i given by (PT2l) 
is exponentially distributed [31] as 

f {it n ) = — exp _- !-!^ , (14) 
7o V 7o / 

and of n is the n th diagonal entry of R7 1 , i.e. 

°i, n = 4^ le M, , (15) 
where e n ^ is the column of I d . . Next we bound of n using Lemma @] 

Lemma 4. Assume Rt w a Hermitian PSD matrix, F is a random N t x d matrix such that 
F*F = Li and R = E {F*R t F}. Then, a 2 n = e^ d R _1 e njrf can be bounded as 

£ < o 2 n < i d = 1 

i (At-A.y)' <ff2< i h t | A, | ,\ , Vne{l,... ,d}, (16) 

Ai "T" AiCAiAjv-dHRtl^) - a n - 4Ai ^ Ajv Ai "T Z J " ^ 1 

where Ai < A 2 < • • • < A at are f/ze eigenvalues o/R t . 

Proof: The Kantrovich inequality 11331 states that for a d X d PSD matrix R [|34l 

^<- W7 ^-(^ + ^ + 2), (17) 
4R(n, n) A ^ / 

where A/6 and A„b are lower and upper bounds for the smallest and largest eigenvalue of R 
respectively. Moreover, R(n, n) = E {F*(n, :)R t F(:, n)} is a quadratic form and can be bounded 
as 

Xi <R{n,n) < X N , (18) 

where Ai < A2 < . . . < \n are the eigenvalues of H t . In addition, since F*F = Lj, using 
Pioncare's separation theorem ll29l . the sorted eigenvalues of R, Aj i G {l,...,d}, can be 
bounded as 

A* < Aj < \ N -d+i i = 1, • • • , d . (19) 
Substituting the bounds from (fl~8l and (fT9l into (fT71) gives an upper bound on o 2 , i.e. 

Note that when d = 1, (fT3l will simplify and (TTSl gives a tighter upper-bound on o\. Combining 
the bounds from (fl~8l and (l20l) results in 



o^ ^ Al ^ rf = 1 > Vne{l,...,d}. 



s 



As Y^ m =i II R( n ) m ) III < ^ II 111 Il35l ' by using the upper bound on the diagonal entries 
of R given by (fTST ), a lower bound on o\ can be found 11361 



\2 



°"n — T~ T77T~~7Tp ^ VnG{l,...,d}. (21) 



Ax Ai(AiAtv — cf R( 



When d = 1, as || R t |||= A^, the right hand side of (12TT) simplifies to 

1 | (Ai-A^y) 2 _ \ N (\ 1 -\ N ) + (\ 1 -\ N ) 2 _ 1 
Ai AiAat(Ai — \n) AiAat(Ai — Ajv) Aat 



which agrees with (|T81) . ■ 
Instead of bounding (fT5l) . we could directly compute Rj. To compute Rj, correlation (covari- 
ance) function between the elements of Fj is needed. In the alternating minimization algorithm, 
the columns of Fj are set to the d { least dominant eigenvectors of 

^ Hfci^AV - C k Cl)H ki , (22) 

k^i 

and because is not independent of H ki for i ^ k, (1221) is not a Wishart matrix and the known 
results for the covariance matrix of the eigenvectors of Wishart matrices (e.g. see 11371 . Il38l and 
references therein) do not lead to an exact characterization of Rj. Assumption \T\ however, paves 
the way for directly computing (approximating) Rj using Corollary Q] and Lemma [51 

Assumption 1. Any set of basis of the interference subspace at the k th receiver, columns of C k , 
is independent of the channels between the interfering transmitters and the k th receiver, Hki far 

i^ke{l,...,K}. 

Using Assumption [T] and Lemma [3l each term of the summation in d22l) is an iV r x iV r Wishart 
distributed matrix with d k degrees of freedom and covariance matrix H t . Moreover, as H^j is 
independent of H nm for k ^ i or i ^ m, based on [30, Theorem 3.4.3], distribution of (|22~1) is 
given by Corollary [Q 

Corollary 1. In the alternating minimization algorithm, using Assumption [7] Lemma [2] and 
Lemma \3\ ft22\) is an N r x iV r complex Wishart matrix with Ylk^i degrees of freedom and 
covariance matrix Hf 

Based on Corollary [H columns of Fj are a subset of eigenvectors of a matrix with central 
complex Wishart distribution. Lemma \5\ gives the covariance (correlation) function between 
components of the eigenvectors of a Wishart matrix which can be used to compute (TT3l . 

Lemma 5. Let X p and u p , for p = 1, . . . , N, be the eigenvalues and the eigenvectors ofanNxN 
matrix with a central complex Wishart distribution with D degrees of freedom and covariance 
matrix H t . Next, let \ p and u p , for p = 1, . . . , N, be the eigenvalues and eigenvectors ofH t - Also 
let the q th elements ofu p and u p be u pq and u pq respectively. Assuming Ai > A 2 > • • ■ > A^ > 0, 
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for D — > oo and fixed N 

E{u p } = u p , (23) 



f \ N \ 



COV(U pq ,U p/q/ J 



n 7a — 3T^ p = p 

D -l,r^ ^ ~ Ap} . (24) 



\p\p/UpqU, p/q/ 



D(X P — A p /) 2 

Proof: See j39l _Chap. 3.6], (371 Chap. 10] and references therein. ■ 
Using Lemma [51 Rj(n, m) Vn, m G {1, . . . , di} in (TTSl is given by 

Ri(n,m) = tr(E{F*(n,:)R t F i (:,m)}) = tr (E {Fi(:, m)F*(n, :)}R t ) 

= tr ((cov (Fi(:, m), F ; (:, n)) + E {Fi(:, m)} E {F*(n, :)}) R t ) . (25) 

Moreover, values given by (1251) should be normalized by noting that F*Fj = 1^ and hence 

tr (cov (Fi(:, m), F^:, m)) + E {F ; (:, m)} E {F*(m, :)}) = 1 Vm G {1, . . . , di}. 

Note that the asymptotic expressions presented in Lemma [5] are valid when the number of 
streams in the IA network is asymptotically large while the number of antennas in each node 
is kept constant. We know, however, for a K-user interference channel, when the number of 
streams in the network increases, IA remains feasible only if the number of antennas at each 
node increases accordingly lfT5l . But such asymptotic results fall into the category of asymptotic 
behavior of eigenvectors of large covariance matrices (e.g. see ll38l and references therein) where 
closed form solutions similar to (|24l) are either not available or their complexity will not benefit 
the current discussion of this paper. It should be noted that more accurate results for the values 
of Rj can always be obtained by using a better approximation or a more complex closed-form 
expression. 

Under certain models of spatial correlation, for example those based on scattering clusters 
[40], the spatial correlation is a function of cluster locations and will vary for different transmit 
and receiver locations. Now we assume there is a potentially different transmit correlation for 
each link pair, i.e. K tz ^ Rt fe for k ^ i G {1, ... , K}. Using Assumption [Q (|22l) will be a linear 
sum of Wishart matrices with equal degrees of freedom but unequal covariance matrices. It is 
shown in iHTI that such a linear sum can be approximated with another Wishart matrix whose 
degrees of freedom and covariance matrix are given by Lemma [6j 

Lemma 6. Suppose T = ^f^ 1 T» where Tj is a N x iV Wishart distributed matrix with di 
degrees of freedom and covariance matrix R ti . Then T has, approximately, Wishart distribution 

■ili rr i i ■ fS 7 J tr(Efi"I 1 d i Rt i ) 2 +tr 2 te^-"I :l diRt i ) 

with a degrees of freedom and covariance matrix R t where a = „ A --i — 



and R t = 3X^=1 ^R ti . 



Proof: See gfl Section 3] ■ 
Using Lemma \5\ and Lemma [6l the correlation between the elements of eigenvectors of d22l) 
can be found and consequently one can compute (fT3l) for each receiver. Note that if H tt 7^ Rt fe 
for i k G {1,...,K}, then af n ^ a\ n . For the rest of this manuscript, to simplify the 
notation, we assume equal transmit correlations at each transmitter. The forthcoming analysis 
and the derived equations, however, can be generalized to unequal transmit correlation matrices. 
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IV. Imperfect Channel Knowledge 

The channel state is estimated and thus known imperfectly in realistic wireless systems. We 
model imperfect CSI through a Gauss-Markov uncertainty of the form 021, ll43l 



H w = y/l - /3m w + /3E , (26) 

where LP ~ CAf(Q, I) is the true Gaussian part of the channel matrix, LP ~ CAf(0, 1) is the 
imperfect observation of H 1 " available to the nodes and E ~ CJ\f(0, 1) is an i.i.d Gaussian noise 
term. The parameter j3 characterizes the partial CSI since (3 = corresponds to perfect channel 
knowledge and (3 = 1 corresponds to no CSI knowledge and values of < (3 < 1 account for 
partial CSI. Note that our forthcoming discussion on the general expression of the imperfect 
CSI at the nodes as modeled in (l26l) can be used to study specific scenarios of imperfect CSI 
such as channel estimation error or analog feedback. The difference between such scenarios is 
how (3 changes as a function of different system parameters. For example, with MMSE channel 
estimation, (3 is a function of pilot symbol SNR PBl or for an analog feedback link, (3 is a 
function of the number of channel uses per channel coefficient and the SNR of the feedback 
link [45 ] . We assume (3 is a constant but one can extend the results to (3 being an arbitrary 
function of the system parameters. 

In presence of a transmit correlation, R t , using arguments similar to [46], IflTll . we assume that 
the transmit correlation varies slower than the channel itself such that the nodes have a perfect 
estimate of H t even if their observation of H is not perfect. Therefore, the channel correlation 
given by (flTT) can be considered in conjunction with the CSI imperfection using the following 
model 

H = (y/l -f3 2 H w + (3E)R] /2 . (27) 

ng and combining matrices are designed i 
effectively ignoring the introduced CSI imperfection. Similar to CQ> an d © 

(V 1 - £ 2 H& + /3E, fc ) R t 1/2 F fc x fc + z 



1/2 

We assume that the precoding and combining matrices are designed using LPR t , thus 
; in 

K 

E 

k=l 

,1/2 



X 



^2k^i AifcXfc 



+ /3E i (R t 1/z ®I x )Fx + z i , (28) 



where H; = [H^R^ 2 F i; Q], A lk is found by solving dA ik = H^R t 1/2 F fc , F and E< are block 
diagonal matrices with Fi, . . . , F^ and Eji, • • • , E«x on their main diagonals respectively and 
x = [xf , . . . , x^-] T . The ZF equalizer, similar to ©, is given by BjH,^ 1 and the distribution of 
the post-processing SINR after applying the ZF equalizer to (1281) is given by Lemma U\ 

Lemma 7. Given the channel model of KT7\) . assuming the IA precoders are designed ignoring 
the CSI imperfection, the post-processing SINR distribution after applying a ZF equalizer is 
given by 



f ci (%, n ) = l/a i exp(- 7iiri /a i ) i e {1, . . . , K}, n e {1, . . . , dj, (29) 



where an = ? ^S x +± y of, n is 8 iven b y (E3 and % = Yn=i tr(F*R t Fi)/di. 



11 



Proof: Applying the ZF equalizer to (1281) , the per stream SINR of the i th user is given by 



P 2 )/d t 



B, (pm^E t (Rl /2 ® I^)FDF*(R, 1/2 ® I i ,)*E*(Hr 1 )* + ^(HJH,)" 1 ) B* 



(30) 



where D = E {xx*}. As Ej has i.i.d unit-variance zero-mean terms, we have E {AEjBE*C} = 
tr(B)AC for any matrices A,B and C independent of Ej lf3~5l Section 21.2]. Therefore, (1301) 
simplifies to 

(1-/3 2 )M 



B t (/3 2 tr((R t 1/2 ® I K )FDF*(R t 1/2 ® I^Hr 1 ^)* + ^(Hj-Hj" 1 ) B* 

(1-/3 2 )M 



B l (/3 2 tr((R t ®I K )FDF*)(HrH i ) 1 + i(HTH i ) ') BJ 



(3D 



where the equality tr(AB) = tr(BA) was used in (f3TT) . Exploiting the block diagonal structure 
of F and (R 4 ® Ir-) results in 



/ 1,71 



(i - n/di 



B, ((/3>T+i)(B?H i )" 1 ) B* 
Similar to ©, (|32b further simplifies to 

(1 - P)/di 



(32) 



7 Cl 



+ l/ 7o ) 



A,^ - QC*)A- 



(33) 



l/2\*, 



H^)*. The SINR distribution follows from comparing 03 to (fl2l . 



where A< = F*(R, 

For R t ^ I, X can either be approximated using Lemma [5] or bounded using (fT6l) by summing 
over the bounds of diagonal entries of Rj. Note that for H t = I, the values of X and of n in 
(|32|) and (TT3T) will be exact and equal to K and 1 respectively. Moreover, as expected, for (3 = 
and of n = 1, d29l) reduces to CUl- Comparing CEl) to (l29l> , the mean SINR at each stream has 
reduced by a factor of <T 4 2 n (7 /3 2 X + 1)/(1 — (5 2 ) which results in the mean post-processing SINR 
reaching the maximum value of (1 — (5 2 ) / '(of n dif3 2 Z) as j — >■ oo, implying a symbol-error-rate 
(SER) floor and a sum rate cap given /5 ^ 0. Note that when (3 = 0, the mean SINR still scales 
linearly with increasing 7 G and the effect of antenna correlation can be seen as a shift in SER 
or sum rate curves. Moreover, if (3 decreases with increasing j a , provided that H t is not rank 
deficient, there will be no SER floor or sum rate cap and the full multiplexing gain promised 
by I A will be attainable. 



V. Comparison with Point-to-Point MIMO 

In an interference channel, instead of utilizing IA, nodes can access the network resources in an 
orthogonal fashion (TDMA/FDMA). In the resulting parallel point-to-point links, nodes can apply 
any of the traditional MIMO transmission strategies. In this section, we discuss the two most 
common strategies, beamforming and SM, and present the post-processing SNR distributions in 
each case. The post-processing SNR distributions can be used to compare IA to beamforming 
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and SM for a wide range of system performance measures. For the point-to-point MIMO links, 
we assume existence of transmit correlation and imperfect CSI similar to (|27l) . 

A transmitter with CSI can use beamforming to send a single stream to its receiver. Assume 
the transmitter ignores the channel imperfection in (T27l) and treats H W R^ 2 as the true channel. To 
maximize the received SNR, the precoding and combining vectors are set to V(:, 1) and U(:, 1)*, 
respectively, where U and V are given by the singular value decomposition of FPR^ 2 = USV* 

1/2 

and £ is diagonal matrix holding the singular values of H W R/ in descending order. The post- 
processing SNR is then 
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BF 



(l-^)A! 
P 2 V + l/la 



(34) 



where v = E{tr (V*(:, l)R t V(:, 1))} and Ai is the largest eigenvalue of U w B} t /2 (R 1 / 2 )* (U w )* 
with correlated central complex Wishart distribution. Note that columns of V are also the 



eigenvectors of (R^ 2 )* (H tu )*H 1 "R£ / z which has central complex Wishart distribution with 
covariance matrix R 4 and iV degrees of freedom and similar to (TT3T) . Lemma [5] can be used 
to compute v. Therefore, the distribution of the beamforming SNR is given by 

/3 2 )/ c -(R 4 , 7 ) 



,1/2 



r (7) 



i 



(35) 



/3V+l/7o 

where f cw (R t ,j) is the distribution of the largest eigenvalue of a correlated central complex 
Wishart matrix with covariance matrix R t ||48l Section IV]. 

In absence of transmit CSI, SM can be used to convey N data streams to the receiver. Assume 
the channel imperfection is ignored at the receivers and the ZF equalizer is given by (H^Rj' 
The post-processing SINR per stream can be written as 



.V2\-l 



SM 
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which is distributed as 



(3 2 )/N 



(/3 2 tr(R t )/N+^ 



(Hi 



1/2, 



,1/2 



In \ln) 



1 

— exp 

CO 



-In 



(36) 



where u 



and cr^(n) is the n th diagonal entry of R t 1 . 
63, the performance of beamforming, SM and IA can be compared 



^(n)(/3 2 tr(R t ) 

Using (f32l), ([34]) and 

for a wide range of system performance metrics. One such metric is the achievable throughput 
given by 



Rsumilo, /3,a) = y2 / lo g2(! + 7)/(7) rf 7, 



(37) 



where d is the number of streams in the network equal to J2 i di, N and 1 for IA, SM and 
beamforming, respectively. Another point of comparison between IA and point-to-point MIMO 
links is the per-stream SER. For any modulation Ai, the AWGN SER is a function of 7, e.g. 
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£ M (j), and the per-stream SER can be written as 



SER A4 (7 ,/3,a) = / £^(7)1(7^7, 



Jo 



where / is any of the distributions given by d32l> . (|34l) or d36]). Note that both ([32]) and ([36]) are 
exponentially distributed and the mean value of the distributions suffices to compare IA and SM 
in terms of per-stream SER (assuming di = dj Vz, j E {1, . . . , K}). Specifically, the ratio of the 
mean SINR of the SM link to the mean SINR of the IA network is given by 



where of n — cr] n for i,j E {1, . . . , K}. When (1381) results in a value greater than 1, the (per 
stream) mean SINR of the SM MIMO network is higher, i.e. when (1381) is larger than 1, given 
a SER constraint, a point-to-point SM MIMO link will satisfy that constraint with a smaller 
transmit power. 



Consider a i^-user N x N MIMO interference channel. All the channel and noise coefficients 
are distributed as CJ\f(0, 1). To focus on the distribution of the received SINR, we exclude 
any large-scale fading effects from the channel model effectively assuming all the nodes are 
co-located and equipped with omni-directional antennas. It is shown in [|9l that using IA, each 
transmitter can send a single stream to its corresponding receiver as long as N = + K)/2\. 
In the forthcoming discussions, the channel coefficient matrices are normalized such that | 
H |||,= N r N t which given the transmit correlation and channel imperfection model of (|26l ) 
translates to tr (R t ) = N t . Note that there exist constraints on the off-diagonal elements of the 
Ri due to the transmit correlation matrix being Hermitian PSD. In this section, for the transmit 
correlation matrix, we adapt a N x N transmit correlation matrix of the form H t (i,j) = a^~^ 
for i,j E {1, • • • ,N}, where a E C is such that H t is positive definite. This model, widely 
used in literature and industry ||49ll , models the correlation between elements of a uniform linear 
array antenna where a = and \a\ = 1 correspond to no correlation and a rank 1 channel, 
respectively. 

The numerical values of R^ 1 computed from (fT3~1) , the theoretical approximation found by 
replacing (1241) and (1231) into (fl"3l and the bounds given by (fT6l) versus a for two different IA 
networks are depicted in Fig. [2l The two IA networks are a 3-user 2x2 and a 5-user 3x3 MIMO 
networks with — 1, i E {1, . . . , K}. Both networks use the alternating minimization algorithm 
to compute the precoding and combining matrices. As expected, the theoretical approximation is a 
better estimate of the true values of R~ l than the presented upper and lower bounds, especially 
for small values of a. Note that by increasing iV and a, the accuracy of the approximation 
decreases. 

In computing (|29l) , for a fixed j a , as (3 increases, the numerator approaches zero and the 
impact of any errors in approximating of n on (|29l) diminishes. Therefore, we expect that by 
increasing the distribution of the received SINR given by (|29l) will approach to the true 
distribution. Moreover, for a fixed (3, as l/j approaches zero, any errors in computing of n will 
be attenuated by a smaller value and we expect (|29l ) to become closer to the true value of the 
distribution. Although there exists numerous methods of quantifying the accuracy of (|29l) , we 



E { f SM ] 




(38) 



VI. Numerical Results 
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choose the Kullback-Leibler divergence (KLD) [ 50] between the distribution given by (1291 and 
the empirical one. The continuous distributions required to compute KLD are approximated with 
uniformly spaced samples in their range. 

Now consider a 3-user IA network with di — 1 Vi 6 {1, 2, 3} and transmit nodes equipped with 
a uniform linear array of 2 antennas spaced a\ apart where a G IR + and A is the transmission 
signal wavelength. Here, we control a by varying a. We use the method proposed in [49] 
(Suburban Macro-cell environment) to calculate a for various antenna spacings. Table IVIIl shows 
the selected antenna spacings and the corresponding a values. The KLD between empirical 
distribution of the received SINR and theoretical distribution given by ([291 as a function of a 
for three values of (3 and two values of 7 Q is shown in Fig. [3] As expected, the divergence 
between the distributions decreases by increasing (3 and j . 

As discussed in Section [TV] the distribution of per stream SINR given by ( 1291 is exact for 
a = 0. The sum rate of a 4-user 3x3 MEMO system, when H t = I, for four values of (3 versus 
7o is depicted in Fig. [4] where theoretical sum rate curves are found by replacing (|29l) into (1371 ) 
and the maximum achievable sum rates (when > 0) are found by replacing the mean value of 
([291 with its limit as 7 D approaches infinity, i.e. / ,^ 2 As can be seen, the numerical results 
exactly follow the theoretical predictions. Moreover, when j3 ^ the multiplexing gain is zero. 
For (3 7^ 0, however, by comparing the sum rate curve with the case of perfect CSI ((3 = 0), one 
could find ranges of transmit power where channel imperfection has practically no effect on the 
sum rate; for example, 7 G < 20dB for (3 = 0.01. 

The effect of transmit correlation and imperfect CSI on the sum rates of a 3-user 2x2 MIMO 
LA network and a 2 x 2 point-to-point MIMO beamforming link is shown in Fig. [51 Several 
important conclusions can be drawn from this figure. First, for non-ideal transmit correlation 
and imperfect CSI, sum rate of a beamforming link can be higher than the sum rate of an IA 
network with the cross-over point between the sum rate curves being a function of a, (3 and 7 D . 
Second, the accuracy of (|29l , as discussed before, increases with increasing 7 Q while it decreases 
with increasing a. Note that as a increases, H t approaches a rank-deficient matrix (violating 
our assumptions in Lemma [5]) which reduces the accuracy of (l29l . Third, Lemma [5] can be used 
to accurately predict the performance of beamforming in presence of transmit correlation and 
channel imperfection. 

Consider a 3-user 2x2 MIMO IA network with d { = l,i e {1,...,K} and a 2 x 2 SM 
MIMO link with a ZF receiver. Although it can be shown that the -R sum of the SM link, in this 
configuration, is always less than the -R SU m of the IA network, a measure of relative performance 
between the two networks can be defined as (f38b . The contours curves of (1381 ), in terms of a and 
(3 for varying 7 Q , are shown in Fig. [61 As can be seen, the SM link has a better performance in the 
areas defined by large (3 and small a and this area grows by increasing 7 D . This behavior can be 
explained by noting that the effect of imperfect CSI on the IA and SM networks is proportional 
to K and N, respectively, and when K > N, imperfect CSI is more destructive for the IA 
network. The effect of antenna correlation, however, is more tolerated in the LA network where 
the eigenvalues of R t , unlike the SM MIMO link, are not directly limiting the ZF performance. 
The corresponding numerical curves of Fig. [6] are shown in Fig. [7J where by increasing j a , the 
theoretical approximation better estimates the true behavior of the IA network. Moreover, both 
the theoretical and numerical contours show similar values for a and (3 for which by increasing 
7o, the relative behavior of the SM link to IA network will not change, i.e. a ~ 0.58 and 
(3 ~ 0.04, showing how our derived theoretical results can be used to accurately predict the IA 
system performance over a large range of (3, a and 7 Q . 
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VII. Conclusion 

The performance of MIMO IA networks depend on the accuracy of CSI and channel correla- 
tion. This paper quantified the impact of imperfect CSI and transmit antenna correlation via the 
per-stream post-processing SINR distribution. Upon using zero-forcing equalizers in a Rayleigh 
channel, post-processing SINR was shown to be exponentially distributed with the mean value 
being a function of the number of antennas at each node, the transmit antenna correlation, the 
imperfection in CSI, and the transmit power. It was shown that, in the presence of imperfect CSI, 
the performance of IA degrades with increasing total number of streams in the network and if 
the imperfection does not vanish at asymptotically high transmit powers, the multiplexing gain 
is zero. Moreover, it was shown that as long as the channel matrices are full-rank, the impact of 
transmit correlation is less detrimental confined to a constant power loss which does not decrease 
the multiplexing gain achievable through IA. The performance of the two most commonly used 
transmit techniques in orthogonal access networks, beamforming and spatial multiplexing, was 
compared to the performance of IA by utilizing the derived SINR distributions where it was 
shown that IA is not always the optimum transmission strategy given realistic system parameters. 

The results of this paper can be used as a starting point for further research. The performance 
of other communication techniques for interference channels such as combinations of orthogonal 
access, beamforming, coordinated multipoint transmission/reception and so on can be compared 
to the performance of MIMO IA networks and optimal switching points between different 
methods based on number of nodes, number of antennas, transmission power, amount of CSI 
imperfection or structure of the transmission correlation can be found which might help guide 
which techniques are appropriate for different networks. 
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Fig. 1. A 7f-user MIMO interference channel where transmitters and receivers have Nt and N r antennas respectively. 




Fig. 2. Comparison of the numerical values of with the proposed approximation and the bounds in Section [TTT] in two 
MIMO IA networks, a 3-user 2x2 and a 5-user 3x3, for varying antenna correlation values. As expected, the approximation 
is closer to the true value compared to the bounds. The proposed approximation is within 10% of the true value for a < 0.3. 
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Fig. 3. Kullback-Leibler divergence between the theoretical distribution of the per-stream received SINR given by l |29l > and 
the true distribution found through numerical simulation as a function of a for varying j3 and j . As can be seen, the accuracy 
of J29l > increases with increasing j3 and transmit power. 



45 




0' ' ' ■ ■ ' 1 

5 10 15 20 25 30 

Transmit Power (dB) 



Fig. 4. Theoretical predictions and numerical results for the sum rate of a 4-user 3x3 constant channel MIMO IA (a = 0) 
for 4 values of /3 versus 7 D . As can be seen, numerical simulations closely follow the analytical curves and the theoretical 
maximum sum rates correctly predict the upper bounds for the given ft. 
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Fig. 5. Achievable sum rate of 3-user 2 x 2 IA MIMO network and a 2 x 2 MIMO beamforming link as a function of j for 
varying a and a fixed ft = 0.19. The accuracy of {29} increases with increasing j and decreases with increasing a. As can be 
seen, beamforming can outperform I A for a wide range of parameters. Moreover, Lemma [5] can be used to accurately predict 
the beamforming performance. 




Fig. 6. Theoretical unity level contour plots for the ratio of the per stream mean SINR in a single-user 2x2 SM MIMO link 
to per stream mean SINR of a 3-user 2x2 MIMO IA network given by d38b for varying f3, 7 and a. As expected, the SM 
link has a better performance for large values of /3 and small values of a. 
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Fig. 7. Numerical unity level contour plots for the ratio of the per stream mean SINR in a single-user 2x2 SM MIMO link 
to per stream mean SINR of a 3-user 2x2 MIMO IA network given by ( 1381 ) for varying f3, y and a. As expected, the SM 
link has a better performance for large values of /3 and small values of a. 



TABLE I 

Channel correlation values for varying antenna spacing in a Sub-urban Macro-cell environment ||49l 



Antenna Spacing (xA) 


a 


a 


10 


-0.1743 + 0.0951? 


0.1986 


9 


0.2064 + 0.1066?' 


0.2323 


8 


-0.0341 - 0.2872?' 


0.2892 


7 


-0.2817 + 0.2408? 


0.3706 


6 


0.4551 + 0.1317?' 


0.4738 


5 


-0.1717- 0.5660?' 


0.5915 


4 


-0.4616 + 0.5439* 


0.7134 


3 


0.8193 + 0.1101?' 


0.8267 



